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Abstract The oxygen quenching effect in Linear Alkl Benzne (LAB) based liquid scintillator (LAB as the 
solvent, 3 g/L 2, 5 diphe-nyloxazole (PPO) as the fluor and 15 mg/L p-bis-(o-methylstyryl)-benzene (bis-MSB) 
\ as the A-shifter) is studied by measuring the light yield as the function of the nitrogen bubbhng time. It is 

■ shown that the light yield of the fully purged liquid scintillator is increased by 11% at the room temperature 

' and the room atmosphere pressure. A simple nitrogen bubbling model is proposed to describe the relationship 

OA . between the relative light yield (oxygen quenching factor) and the bubbling time. 
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1 Introduction 

It has been shown experimentally and theoreti- 
cally that the presence of oxygen in the liquid scintil- 
lator (LS) can lower the light yield, modify the fluo- 
rescence pulse shape, shorten the attenuation length 
and decrease the positron annihilation lifetime'^'— 
For most aromatic molecules, the quenching of the 
electronically singlet state (Si) leads to the formation 
of triplet state (Ti). The oxygen molecule is some- 
what special. Its ground state is a triplet and the 
the next state is a singlet lying about 0.98 eV over 
the ground state'"''; Oxygen molecules in aromatic 
molecules can absorb the energy of singlet state of 
aromatic molecules, and make the spin allowed tran- 
sition to the triplet state. This decreases the fluo- 
rescence. Such a transition only occurs in aromatic 
molecules which have a energy gap of Si -Ti greater 
than 0.98 eV. Most of aromatic molecules meet this 
requirement. In most experiments, the dissolution 
of oxygen into LS is undesirable. Since this brings 
uncertainties to the experiments. Hence, the oxygen 
quenching effect should be well studied and deter- 
mined. 

Usually, there are three ways to eliminate the dis- 
solved oxygen fromsolutions: 1) the vacuum distilla- 
tion; 2) the ultrasonic^; 3) the nitrogen(or argon) 
bubbling'^-. In neutrino experiments, large quantity 
of liquid scintillator is required. The most economical 



and practicable way to eliminate oxygen in LS is the 
nitrogen (or argon) bubbling. 

Linear alkyl benzne (LAB), which is composed of 
a linear alkyl chain of 10-13 carbon atoms attached 
to a benzene ring , is a low cost product of petro- 
chemical industry and is often used as the material 
of detergent. Its aromatic structure makes it be use- 
ful as a scintillator solvents inherently. It has many 
appealing properties, including the high flash point 
(130 °C), the low toxicity, the high light yield and ex- 
cellent transparency'^. LAB based liquid scintillator 
will serve as the antineutrino target in the Daya Bay 
neutrino experiment'^'. 

In this work, we measured the effect of oxygen 
on the light yield of LAB LS, and built a nitrogen 
bubbling model to describe the relationship between 
relative light yield (the oxygen quenching factor) and 
the nitrogen bubbling time. Parameters in the model 
were determined by our experimental data. 

2 Nitrogen bubbling model 

When LS is exposed to the air, oxygen molecules 
dissolved in the LS exchanges with those in the air. 
This process is in dynamical equilibrium. It is rea- 
sonable to assume that the oxygen dissolved in the 
un-bubbled LS is saturated due to its long time con- 
tact with the air. This means that the number of 
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oxygen molecules dissolving into the LS is equal to 
those escaping from the LS. When LS is flushed with 
nitrogen, the oxygen partial pressure in the nitrogen 
bubble, which presents in LS, can be thought to be 
zero. Therefore, oxygen molecules will escape from 
LS and enter the nitrogen bubbles. The dissolution 
of oxygen molecules into the LS can be ignored in 
the interface of nitrogen and LS. Since that the oxy- 
gen molecules diffusion rate into LS is much higher 
than the oxygen escaping rate, which denotes that 
oxygen molecules are uniformly distributed in LS; 
It is reasonable to assume that the oxygen escaping 
rate is proportional to the contact area of bubbles 
with LS and the oxygen partial pressure in LS which 
is proportional to the oxygen concentration. Then, 
the equation describing the variation of the oxygen 
molecule number in LS, dN/dt, is given by 



dN 
~dt 



(1) 



where [Q] is the oxygen concentration dissolved in LS; 
fee and 5" are the oxygen escaping rate and nitrogen- 
LS contact area, respectively. The oxygen concentra- 
tion can be expressed as the following: 



N 



(2) 



where N and are the oxygen molecule number in 
LS and the LS volume, respectively. Hence, Eq. ([T]) 
can be rewritten as 



dN 



it' 



(3) 



Then, the variation of the oxygen molecule num- 
ber in LS is 

dN KS , 

= —dt. (4 

Eq. can be written in the integration form: 

k„S . 



7V = 7Voexp(-^t), 



(5) 



where A'o is the oxygen molecule number in LS with- 
out bubbling. 

At low concentration, the quenching of fluores- 
cence by a quencher in solution can be described by 
the well-known Stern- Volmer relationship. 



(6) 



where Iq is the intensity or rate of fluorescence with- 
out a quencher present, / is the intensity or the rate 
of fluorescence with a quencher, [Q\ is the quencher 
concentration dissolved in LS, and kq is quenching 
constant. 



Hence, the relative light yield as the function of 
the bubbling time is given by 



koNo 



^^1 



■ exp 



/ KS 



(7) 



Note that Nq/Vs is the saturated concentration of 
oxygen which is dependent of temperature and at- 
mospheric pressure; At fixed temperature and atmo- 
spheric pressure, iVo/K is a constant. Denote the 
kqNo/Vs item by a constant A. Then, Eq. ([7]) turns 
to be 

^0 , . < / k^S , 



(8) 



The oxygen quenching factor, fq, is defined as the 
light yield of LS with oxygen to that without oxygen, 
i.e. I/Iq. Eq. ([5]) can be expressed in the form of /q. 



/p = l/(l + Acxp(-Mt)). 



Nitrogen 



(9) 
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Glass tube 



Fig. 1. Bubbling setup 

The bubbling setup can be illustrated by the sim- 
plified plot in Fig. [T] The nitrogen-LS contact area 
consists of two parts: the area of the nitrogen bubble 
surface and the contact surface area of liquid level. 



S ~ nSt + Si. 



(10) 



where n is the average number of nitrogen bubble 
present in the LS. 5*^ and Si are the average surface 
area of bubbles and the area of the liquid level, re- 
spectively. Then, Eq. ^ can be rewritten as 



/q = 1/(1 + AeM~k,{nSt + Si)t/V,)). 



(11) 



The parameters A and fc^ are constants depend- 
ing on the temperature and pressure for specific LS. 
In the following sections we will evaluate A and k^ 
for LAB LS at the room temperature and the atmo- 
spheric pressure experimentally. 
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3 Experiment 

In order to observe the light output variation due 
to the oxygen quenching, sets of samples of LS (LAB 
as the solvent, 3 g/L PPO as the fluor and 15 mg/L 
bis-MSB as the wavelength shifter) were bubbled with 
different time. Fig. [T] shows the bubbhng setup. Six 
sets of LS samples (50 ml for each sample) were used. 
One set of samples was not bubbled, and the other 
five were bubbled with high purity nitrogen for 5 min 
(200 ml), 12.5 min (500 ml), 18.7 min (750 ml), 25 
min (1000 ml) and 31.25 min (1250 ml), respectively 
(the numbers in the brackets arc the nitrogen vol- 
umes). Nitrogen flow rate was precisely controlled at 
40 ml/min by a flowmeter. The bubbles present in LS 
can be thought to be spherical. The bubble number 
appears in the tube can be easily counted and the 
diameters of bubbles can be measured with rulers. 
In our condition, 4 bubbles, with 4 mm in diameter, 
present in LS. The diameter of the tube is 23 mm. 
Hence, the LS-nitrogen contact area is 616.5 mm^. 



Fig. 2. Schematic view of the experimental setup. 

The setup shown in Fig. [2] was used to measure 
the light output. LS was encapsulated in a cylindrical 
teflon cell (5 cm in diameter and 2.5 cm in height). 
The end of the cell was terminated with UV glass 
which was coupled to a 2-inch high energy resolution 
PMT (Hamamatsu CR105). The cell and PMT were 
placed in a dark steel box. The cell was exposed to a 
^^''Cs 7 -ray source. 



Fig. 3. System stability in 12 hours. The sta- 
bility was tested by a LED driven by a pulse 
generator 



The charge of PMT output pulse was firstly cov- 
ersed to amplitude by a charge sensitive amplifier and 
then shaped by shaping filter. Finally, the shaped 
signal amplitude was analyzed by a pulse height ana- 
lyzer. System stability has been measured by means 
of a LED driven by pulser generator. Fig.[3]shows the 
system stability. System trended stable after 6 hours 
burning. System was burned about 12 hours before 
data acquisition in our experiment. 

4 Data Analysis and Results 




Deposited Energy (KeV) 

Fig. 4. Energy response of LS to the ^''^Cs 7- 
ray source simulated by GRESP. The counts 
are normalized to the number of 7-rays emit- 
ted from the '^^'^Cs source. 

Energy spectrum N{E) of the Compton scatter- 
ing electrons is generated by means of Monte Carlo 
GRESP code. Fig. [D shows the simulation result. It 
should be noted that the resolution smearing is not 
considered in the GRESP Monte Carlo code. We con- 
sidered the resolution smearing similar to the way de- 
scribed in Ref 0| . The "realistic " Monte Carlo spec- 
trum can be obtained from the convolution of the 
simulation spectrum with the system response func- 
tion. 



\H)^ j R{H,L)NL{L)dL, 



(12) 



where R{H,L) is the response function, H is the ADC 
channel and Ni^{L) is the spectrum of light output. 
For electrons (above about 50 KeV) , the light output 
L, the energy emitted as fluorescence, is proportional 



to the energy E deposited in the LS '. , i.e 



(13) 



where S is the absolute scintillation efficiency. The 
light output L can be defined to be in units such that 
it is equal to the electron energy E'^, i.e 5*= 1 and 



L{E)^E. 



(14) 
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When oxygen presents in LS, the light yield will de- 
crease. Let the quenching factor, i.e. the light yield 
with the oxygen quenching to that without the oxy- 
gen quenching, ///qj be fg. Then the light output 
can be rewritten as 



(15) 



It should be reminded that /q is 1 for LS without 
a quencher. Considering that the spectrum of de- 
posited energy is N{E), the spectrum L, Nl{L), is 
N{L/ fgS)/ fg. Assuming that the response R{H,L) 
of the detector for the fixed light output is Gaussian, 



R{H,L)^Bcxp 



{H-cLf 



(16) 



where c is the light output to ADC channel conver- 
sion factor, and B is the normalization factor between 
the "realistic" Monte Carlo spectrum and the exper- 
imental spectrum. a^L can be expressed in the form 
of the system resolution, 



cL 



(17) 



where p = A{cL)/cL is the detector resolution, i.e. 
FWHM. and the background can be thought to be 
expositional'®' , 



iV^«(i/)=Ciexp -C2F + C3 



(18) 



where Ci, C2, C3 are parameters of the background. 
Then, the expected experimental spectrum can be 
written as 

N'^'{H) = B f exp(- ^^~''{'^^^ )N{E)dE 



with 



Oe = 



-hciexp^caij-l-ca), (19) 
c/qEp 



(20) 



The expected experimental spectrum is determined 
by _B,Ci,C2,C3, c, /q and the detector resolution p. c 
depends on the optical properties of the cell, PMT 
and electronics. It is independent of quenching and 
can be treated as the a constant in our experiment. 
Hence the item fqc, which is proportional to the 
quenching factor, can be taken as one free parameter. 
To evaluate the free parameters we used the ROOT 
package (CERN data analysis package) to fit the ex- 
perimental spectrum with Eq. (jl9|) . Fig.[5]shows the 
fit result of the background, the realistic spectrum, 
and the expected spectrum for the experimental spec- 
trum. The items /gc for different LS samples can be 
obtained by fitting their experimental spectrums with 
Eq. (fTO)) . With the increase of the nitrogen bubbling 
time, the light output increased. We assume that 
oxygen was fully bubbled, namely no oxygen quench- 
ing effect remains in LS, when the light yield change 
little. It should be reminded that the quenching fac- 
tor, ///o, for full bubbled LS is 1. c was established 
from the fit result of the full bubbled LS. Then, the 
value of /q for different LS samples were determined. 
Fig. [6] shows the /g for six LS samples. The solid 
line in Fig. [H] shows the fits with Eq. . From the 
experimental result, we know that the light yield is 
increased by about 11% (20°C) by means of removing 
the oxygen in LS. 



X' / Mlf 0.9725 / 4 

A 0.I284 + 0.01I81 
K 1.134 ±0.07997 
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Fig. 5. (1) Experimental spectrum, (2) "real- 
istic" Monte Carlo spectrum ,(3) Background 
defined in Eq. (|18[l. (4) expected experimental 
spectrum defined in Eq. (fTU|l for ^^'^Cs. 



Fig. 6. /g as a function of the bubbling time. 
The solid curve shows the fits with Eq. ((9)1. 



5 Conclusions 

The oxygen quenching in LAB liquid scintillator 
and the degassing model have been studied. From 
the experiment, we know that LAB LS light yield is 
increased by 11% by fully removing oxygen at 20° C. 
Moreover, we proposed a model to determined the 
relationship between the light yield and the bubbling 
time in this paper. The parameters in the model have 
been fixed experimentally. 
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